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The identification of computational and experimental reduced-order models 
(ROMs) for the analysis of unsteady aerodynamic responses and for efficient 
aeroelastic analyses is presented. For the identification of a computational aeroe- 
lastic ROM, the CFL3Dv6.0 computational fluid dynamics (CFD) code is used. 
Flutter results for the AGARD 445.6 Wing and for a Rigid Semispan Model 
(RSM) computed using CFL3Dv6.0 are presented, including discussion of associ- 
ated computational costs. Modal impulse responses of the unsteady aerodynamic 
system are computed using the CFL3Dv6.0 code and transformed into state-space 
form. The unsteady aerodynamic state-space ROM is then combined with a 
state-space model of the structure to create an aeroelastic simulation using the 
MATLAB/SIMULINK environment. The MATLAB/SIMULINK ROM is then 
used to rapidly compute aeroelastic transients, including flutter. The ROM shows 
excellent agreement with the aeroelastic analyses computed using the CFL3Dv6.0 
code directly. For the identification of experimental unsteady pressure ROMs, 
results are presented for two configurations: the RSM and a Benchmark Su- 

percritical Wing (BSCW). Both models were used to acquire unsteady pressure 
data due to pitching oscillations on the Oscillating Turntable (OTT) system at 
the Transonic Dynamics Tunnel (TDT). A deconvolution scheme involving a step 
input in pitch and the resultant step response in pressure, for several pressure 
transducers, is used to identify the unsteady pressure impulse responses. The 
identified impulse responses are then used to predict the pressure response due 
to pitching oscillations at several frequencies. Comparisons with the experimental 
data are presented. 
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Introduction 

A T present, the development of CFD-based 
reduced-order models (ROMs) is an area 
of active research at several industry, govern- 
ment, and academic institutions . 1 Development of 
ROMs based on Volterra theory is one of several 
ROM methods currently under development . 2 ~ 7 
Reduced-order models based on Volterra theory 
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have been applied successfully to Euler and Navier- 
Stokes models of nonlinear unsteady aerodynamic 
and aeroelastic systems. Volterra-based ROMs are 
focused on the creation of linearized and nonlin- 
ear unsteady aerodynamic impulse responses that 
may then then used in a convolution scheme to 
provide the linearized and nonlinear responses of 
the system to arbitrary inputs. In this setting, the 
linearized and nonlinear impulse responses are the 
ROMs of the particular nonlinear system under in- 
vestigation. Upon transformation of the linearized 
and nonlinear impulse responses into state-space 
form, the state-space models generated can also 
be considered ROMs. 

Another ROM method, different from the 
Volterra-based ROM approach, is the Proper Or- 
thogonal Decomposition (POD) technique. The 
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POD technique is a method that is used exten- 
sively at several research organizations for the de- 
velopment of reduced-order models. A thorough 
review of POD research activities can be found in 
Beran and Silva. 1 In addition, a review of the is- 
sues involved in the development of reduced-order 
models for fluid-structure interaction problems is 
provided by Dowell and Hall. 8 A topic of re- 
cent interest is the potential development of hybrid 
POD/Volterra methods. These hybrid techniques 
would combine the spatial resolution possible with 
POD methods with the low dimensionality and 
computational efficiency of Volterra methods. 9 

The linearization of a nonlinear aeroelastic 
model is an important first step towards under- 
standing the nature and magnitude of nonlinear 
aeroelastic phenomena. The response of a lin- 
earized system about a nonlinear steady-state con- 
dition can be obtained via several methods. Some 
of these methods include the order reduction of 
state-space models using various techniques . 10, 11 
One method for building a linearized, low-order, 
frequency-domain model from CFD analysis is to 
apply the exponential (Gaussian) pulse input. 12 
This method is used to excite an aeroelastic sys- 
tem, one mode at a time, using a smoothly- 
varying, small-amplitude Gaussian pulse. The 
time-domain aeroelastic responses due to the expo- 
nential pulse input are transformed into frequency- 
domain generalized aerodynamic forces (GAFs). 
These linearized GAFs can then be used in stan- 
dard linear frequency-domain aeroelastic analy- 
ses. 13 Raveh et al 14 applied this method but re- 
placed the exponential pulse input with step and 
impulse inputs. Raveh 15 also performed para- 
metric variations in order to better understand 
the numerical issues associated with impulse and 
step responses, particularly for nonlinear prob- 
lems. Guendel and Cesnik 16 applied the Aerody- 
namic Impulse Response (AIR) technique, based 
on Volterra theory, to the PMARC aerodynamic 
panel code. The PMARC/AIR code was applied 
to a simplified High Altitude Long Endurance 
(HALE) aircraft for rapid linear and nonlinear 
aeroelastic analysis of the vehicle. 

As mentioned above, various inputs can be used 
in the time domain (using a CFD code) to gen- 
erate GAFs in the frequency domain in order 
to perform standard, frequency-domain aeroelas- 
tic analyses. But if time-domain aeroservoelastic 
(ASE) analyses are desired, the frequency-domain 
GAFs are transformed back into the time domain 
using traditional rational function approximation 
(RFA) techniques. These techniques include, for 
example, the well-known Roger’s approximation 17 
and the Minimum State technique. 18 The RFA 
techniques transform frequency-domain GAFs into 
state-space (time domain) models amenable for use 


with modern control theory and optimization. The 
process just described transforms time-domain in- 
formation (CFD results) into frequency-domain 
information only to have the frequency-domain in- 
formation transformed back into the time domain. 

Gupta et al 19 and Cowan et al 20,21 applied a 
set of flight testing inputs to an unsteady CFD 
code and used the information to create a lin- 
ear ARMA (autoregressive moving average) model 
that was transformed into state-space form. Al- 
though this technique is applied entirely within 
the time domain, the shape of the inputs applied 
to the CFD code requires tailoring in order to 
excite a specific frequency range, resulting in an 
iterative process. In a similar vein, Rodrigues 22 
developed a state-space model for an airfoil in tran- 
sonic flow using a transonic small-disturbance algo- 
rithm. In the present paper, a direct approach for 
efficiently generating linearized unsteady aerody- 
namic state-space models is presented. 23 Although 
the present application deals with linearized re- 
sponses based on linearized impulse responses (lin- 
earized Volterra kernels), the method can be for- 
mally extended to address nonlinear aeroelastic 
phenomena via the use of nonlinear impulse re- 
sponses (nonlinear Volterra kernels). 

The development of ROMs for computational 
analyses has proven beneficial to understanding 
dominant flow physics. 1,23 These ROMs provide 
insight regarding the level of nonlinearity within 
a physical process including nonlinear aeroelastic 
responses and aeroelastic limit cycle oscillations 
(LCO). Likewise, the development of ROMs for the 
analysis of experimental data can provide valuable 
insights regarding the nature of the experimental 
flow physics and structural response. 

Experimental investigation of complex flight dy- 
namic and aeroelastic phenomena are best under- 
stood by studying the underlying unsteady aero- 
dynamics. To this end, experiments designed to 
measure the unsteady aerodynamic response of 
various configurations provide significant and valu- 
able information. 24-27 Experimental results are 
compared to various types of numerical analyses 
(such as CFD) to provide insight into the underly- 
ing physics of the problem. 

Insight gained from unsteady aerodynamic ex- 
periments and analyses can then be used to alter 
or control some aspect of the underlying physics 
leading to a modification of the performance of a 
vehicle. This is the primary goal of the flow control 
research effort. 28-31 A thorough understanding of 
the dominant flow physics can lead to an optimal 
flow control strategy. Therefore, understanding 
unsteady aerodynamic behavior is an essential step 
in the design of flow control concepts. 

This paper has two purposes: (1) to present the 
methodology and sample applications of identify- 
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ing computational reduced-order models for aeroe- 
lastic analyses; and (2) to present the methodol- 
ogy and sample applications of identifying exper- 
imental reduced-order models for unsteady aero- 
dynamic analyses. As part of the identification 
of computational ROMs, the paper will discuss 
the development of linearized, unsteady aerody- 
namic state-space models for aeroelastic analysis 
and flutter prediction using the parallelized, aeroe- 
lastic capability of the CFL3Dv6 code. 23 As such, 
this paper begins with a brief description of the 
CFL3Dv6.0 code. This is followed by an outline 
of the ROM development process, and a descrip- 
tion of the CFD-based impulse response technique. 
The Eigensystem Realization Algorithm (ERA), 32 
which transforms an impulse response into state- 
space form, is then described. Flutter results 
for the AGARD 445.6 Aeroelastic Wing and the 
RSM using the CFL3Dv6.0 code are presented, 
including computational costs. Unsteady aerody- 
namic state-space models are then generated and 
coupled with a structural model within a MAT- 
LAB/SIMULINK 33 environment for rapid calcu- 
lation of aeroelastic responses, including flutter. 
Aeroelastic responses computed directly using the 
CFL3Dv6.0 code are compared with the aeroelas- 
tic responses computed using the CFD-based ROM 
within the MATLAB/SIMULINK environment. 

As part of the identification of experimental 
ROMs, the paper demonstrates the feasibility 
of identifying unsteady aerodynamic impulse re- 
sponses from experimental unsteady aerodynamic 
measurements. 34 These experimentally-defined 
aerodynamic impulse responses can be used to pro- 
vide insight regarding the dominant flow physics 
of the experiment. The paper includes a brief de- 
scription of the methods and experimental facilities 
used in the identification of these ROMs. Finally, 
the paper investigates the filtering capability of the 
deconvolution process that is used to generate the 
aerodynamic impulse responses. 

Description of Computational 
Methods 

The following subsections introduce the paral- 
lelized, aeroelastic version of the CFL3Dv6.0 code 
and an outline of the ROM development process. 
Two phases of the ROM development process are 
then described: the identification of unsteady aero- 
dynamic impulse responses and the transformation 
of these impulse responses into state-space form. 

CFL3Dv6.0 Code 

The aeroelastic computational fluid dynamics 
code used in this study is CFL3Dv6.0 code, which 
solves the three-dimensional, thin-layer, Reynolds 
averaged Navier-Stokes equations with an upwind 
finite volume formulation . 35,36 The code uses 
third-order upwind-biased spatial differencing for 


the inviscid terms with flux limiting in the pres- 
ence of shocks. Either flux-difference splitting 
or flux- vector splitting is available. The flux- 
difference splitting method of Roe 37 is employed 
in the present computations to obtain fluxes at cell 
faces. There are two types of time discretization 
available in the code. The first-order backward 
time differencing is used for steady calculations 
while the second-order backward time differenc- 
ing with subiterations is used for static and dy- 
namic aeroelastic calculations. Furthermore, grid 
sequencing for steady state and multigrid and local 
pseudo-time stepping for time marching solutions 
are employed. In this paper, multiblock MPI par- 
allel aeroelastic computations, including flutter, 
are performed. 

ROM Development Process 

An outline of the ROM development process is 
as follows: 

1. Implementation of impulse response tech- 
nique into aeroelastic CFD code; 

2. Computation of impulse responses for each 
mode of an aeroelastic system using the aeroelastic 
CFD code; 

3. Impulse responses generated in Step 2 are 
input into the ERA code; 

4. Evaluation/ validation of the state-space mod- 
els generated in Step 3; 

Steps 1 and 2 are described in greater de- 
tail in the references that address Volterra-based 
Reduced-Order Models (ROMs) such as Refs. 3-7. 
The basic premise of Volterra-based ROMs is the 
extraction of linear and nonlinear kernel functions 
that capture the input-output functional relation- 
ship between, for example, unsteady motion of a 
wing (input) and the resultant loads created by 
that motion (output). For Volterra-based ROMs, 
these kernel functions are linearized/nonlinear im- 
pulse response functions. A brief overview of Steps 
2 and 3 is presented in the following sections. 23 

CFD-Based Discrete Unit Impulse Response 
Technique 

An aeroelastic system can be viewed as the cou- 
pling of an unsteady aerodynamic system (flow 
solver) with a structural system. The majority of 
computational aeroelastic methods can be repre- 
sented by the coupling of the flow solver (nonlin- 
ear aerodynamic system) with a linear structural 
model (Figure 1). The present study focuses 
on the development of an unsteady aerodynamic 
ROM (Figure 2) that is then coupled to a struc- 
tural model for aeroelastic analyses. 

A standard technique for computing linearized 
generalized aerodynamic forces (GAFs) for an 
aeroelastic system with n modes using a CFD code 
is the application of a Greens function (influence 
function) approach. Using the CFD code, each 
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Fig. 1 Coupling of structure and aerodynamics 
within an aeroelastic CFD code. 
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Fig. 2 Identification of generalized aerody- 
namic forces (GAFs). 

structural mode is individually excited to obtain 
the response of all the modes: to this excitation. 
This process is applied to all n modes, resulting in 
an n by n ’’matrix” of responses. The term ’’ma- 
trix” is in quotps to indicate that the responses ob- 
tained using this method are usually time-domain 
functions rather than constants that usually pop- 
ulate a standard matrix. 

This technique is a linearization by virtue of the 
fact, that, in a computational aeroelastic analy- 
sis, the input to the nonlinear flow solver is the 
total physical deformation of the wing consisting 
of the weighted summed total of all the modal 
responses of interest. By applying a separate ex- 
citation to each mode through the nonlinear flow 
solver, the total nonlinear aeroelastic response is 
being approximated by a linear superposition of 
its individual responses. For a linear flow solver, 
this approach would be exact.. This approxima- 
tion is valid only for small input, amplitudes. This 


is not. necessarily a. drawback as, quite often, the 
linearized dynamic aeroelastic response about, a. 
nonlinear steady (or static aeroelastic) condition is 
a. reasonable representation of the nonlinear a.eroe- 
la.st.ic system under investigation. Although not. 
presented in this paper, methods have recently 
been developed that, will enable simultaneous ex- 
citation of all the modes. 

There are three types of modal excitation in- 
puts that, are typically used when implementing 
this technique. The first, is a. brute-force approach 
based on the input, of sine waves of individual fre- 
quencies. The individual modal responses to these 
inputs for n modes and r frequencies requires n 
times r separate code evaluations. In addition, the 
time length required for each one of these eval- 
uations can be quite large (i.e. , computationally 
expensive) in order to get. an adequate number of 
cycles for adequate frequency resolution, especially 
for the lower frequencies. This approach is clearly 
the least, efficient.. 

A second, more elegant, approach, involves the 
use of an Exponential (Gaussian) pulse . 12 The 
exponential pulse can be shaped to excite a. partic- 
ular range of frequencies. Because an exponential 
pulse excites a. pre-selected frequency range, only 
one code evaluation is required per mode. This 
is a. significant, computational savings compared to 
the brute-force approach, but. shape optimization 
of the exponential pulse may be required when 
targeting a. particular frequency range. In addi- 
tion, the exponential pulse appears to be strictly 
limited to linearized analyses. Whereas the im- 
pulse function finds formal application to nonlinear 
problems via. Volt.erra. theory, the inclusion of the 
exponential pulse within a. Volt.erra.-t.ype theoreti- 
cal framework is undefined. 

A third, recently-developed, approach consists of 
replacing the exponential pulse input, with a. unit, 
impulse. 3H:> The unit, impulse excites the entire 
frequency range of a. system so that, shape opti- 
mization to excite a. particular range of frequencies 
is not. necessary. In addition, due to the simplicity 
of the input, and the short, amount, of time required 
for convergence, each solution is computed with 
significant, computational efficiency. 

System /Observer/Controller Identification 
Toolbox (SOCIT) 

In structural dynamics, the realization of 
discrete-time state-space models that, describe the 
modal dynamics of a. structure has been en- 
abled by the development, of algorithms such as 
the Eigensyst.em Realization Algorithm (ERA ) 32 
and the Observer Kalman Identification (OKID ) 38 
Algorithm. These algorithms perform state- 
space realizations by using the Markov parame- 
ters (discret.e-t.ime impulse responses) of the sys- 
tems of interest.. These algorithms have been 
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combined into one package known as the Sys- 
tem/Observer/Controller Identification Toolbox 
(SOCIT) 30 developed at NASA Langley Research 
Center. 

The primary algorithm within the SOCIT group 
of algorithms used for the present system realiza- 
tion is the ERA. A brief summary of the basis of 
this algorithm follows. 

A finite dimensional, discrete-time, linear, time- 
invariant. dynamical system has the state- variable 
equations 


x(k + 1) = Ax(k) + Bu(k) (1) 

y(k) = Cx(k ) + Du(k) (2) 

where x is an n- dimensional state vector, u an 
m- di ma nsion a,1 control input vector, and y a p- 
dimensional output or measurement vector with 
k being the discrete time index. The transition 
matrix, A, characterizes the dynamics of the sys- 
tem. The goal of system realization is to generate 
constant matrices (A, B, C, and D) such that the 
output responses of a given system due to a par- 
ticular set of inputs is reproduced by the discrete- 
time state-space system described above. 

For the system of Eqs. (1) and (2), the time- 
domain values of the discrete-time impulse re- 
sponse are also known as the Markov parameters 
and are defined as 

Y(k) = C'A k ~ 1 B (3) 

with B an (n x m) matrix and C a (p x n) matrix. 
The ERA algorithm begins by defining the general- 
ized Hankel matrix consisting of the discrete-time 
impulse responses for all input/output. combina- 
tions. The algorithm then uses the singular value 
decomposition (SVD) to compute the A, B, and C 
matrices. 

In this fashion, the ERA is applied to unsteady 
aerodynamic impulse responses to construct un- 
steady aerodynamic state-space models. 

Description of Experimental 
Facilities, Models, and Methods 

Transonic Dynamics Tunnel (TDT) 

The Langley Transonic Dynamics Tunnel (TDT) 
is a unique national facility dedicated to identi- 
fying, understanding, and solving relevant aeroe- 
la.stic problems. The TDT is a closed-circuit, 
continuous-flow, variable- pressure, wind tunnel 
with a 16-foot square test section with cropped 
corners. The tunnel uses either air or a heavy gas 
as the test medium and can operate at stagnation 
pressures from near vacuum to atmospheric, has a 
Mach number range from near zero to 1.2 and is 
capable of maximum Reynolds numbers of about 
3 million per foot in air and 10 million per foot 
in heavy gas. The current TDT heavy gas test 
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Fig. 3 Side view of the Oscillating Turntable 
(OTT). 


medium is 1,1, 1,2 t.etra.fluoroet.ha.ne, R-134a. 40,41 
The TDT is specially configured for flutter test- 
ing, with excellent model visibility from the con- 
trol room and a. rapid tunnel shutdown capability 
for model safety (bypass valves). Model mount 
systems include a sidewall turntable for semispan 
models, a variety of stings for full-span models, an 
Oscillating Turntable (OTT), and a cable-mount 
system for ’’flying” models. 

Oscillating Turntable (OTT) 

The OTT is a unique research tool at the TDT 
that provides the ability to oscillate relatively 
large, semispa.n wind-tunnel models in pitch at fre- 
quencies up to 40 Hz. This research tool has been 
designed specifically for the acquisition of unsteady 
pressure and loads data on rigid wind-tunnel mod- 
els in order to study flow phenomena associated 
with flutter, LCO, shock dynamics, and nonlin- 
ear unsteady aerodynamic effects on a wide vari- 
ety of aerospace vehicle configurations at transonic 
speeds. It is anticipated that unsteady pressure 
measurements due to precisely controlled model 
motions will provide valuable data for CFD corre- 
lation and aircraft design with respect to unsteady 
aerodynamic/aeroelastic phenomena.. 25 Models 
may be oscillated sinusoidally at constant, or vary- 
ing frequencies, be subjected to a. step input., or 
undergo user-defined motions. 

Figure 3 highlights key components of the OTT. 
The OTT utilizes a. powerful rotary hydraulic actu- 
ator, rated at 495,000 in-lbf, and a. digital Propor- 
tional, Integral, Derivative, Feedforward (PIDF) 
control system to position and oscillate models. 
Power for the OTT is supplied by a 3000 psi, 150 
gpm hydraulic power unit which is located outside 
the tunnel pressure shell. 

Rigid Semispan Model (RSM) 

The RSM pla.nform is a 1 / 12th scale configura- 
tion based on an early design known as the Ref- 
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Fig. 4 Planform, model details, and instru- 
mentation layout for the RSM wind-tunnel 
model. 

erence H configuration that was a component of 
the High Speed Research (HSR) program. Model 
airfoil shapes were based on those of the Reference 
H, with the model wing thickness being increased 
to a constant 4% thickness-to-chord ratio in order 
to accommodate pressure instrumentation at the 
wing tip. The model was designed to be very stiff 
to allow the measurement of aerodynamic proper- 
ties with only negligible effects of structural defor- 
mations. 

Figure 4 shows the planform layout and main 
components of the RSM including the OTT mount. 
The leading and trailing edges were removable in 
order to access pressure instrumentation in those 
regions. A removable tip cap allowed access to 
pressure instrumentation at the wing tip. The 
RSM could be tested either with or without a pair 
of flow-through nacelles. The nacelles were rigidly 
attached to pylons on the lower, inboard surface 
of the wing. The RSM wing had a graphite epoxy 
composite structure with an open-cell foam core. 

The instrumentation layout for the RSM (visi- 
ble in Figure 4) consisted of 131 insitu unsteady 
pressure transducers located at the 10, 30, 60, and 
95%) span stations. Six additional unsteady pres- 
sure transducers were installed at the 20% chord 
station for the 20, 45, and 75%) span stations for 
both upper and lower surfaces. Channels were 
carved into the foam core to accommodate the 
wiring for the instrumentation. Instrumentation 
also included accelerometers installed throughout 
the wing. The fuselage fairing used for testing the 
RSM on the OTT was instrumented with unsteady 
pressure transducers. For the flutter analysis us- 
ing the RSM, 1 the first four structural modes were 
softened in order to create a flutter response. The 
purpose of this model was to compare full CFD 
aeroelastic analyses with ROM aeroelastic analy- 
ses and to use the data acquired with the RSM on 
the OTT for identification of experimental ROMs. 

Benchmark Supercritical Wing (BSCW) 

The BSCW is shown mounted in the TDT test 
section in Figure 5. The model has a rect.angu- 



Fig. 5 The Benchmark Supercritical Wing 
wind-tunnel model. 


lar planform with a 32 inch span, 16 inch chord, 
and a NASA SC(2)-0414 airfoil. Using 40 in situ 
transducers, unsteady pressure measurements were 
made along the chord at the 60 percent spanwise 
location at Mach numbers ranging from 0.4 to 0.85 
and dynamic pressures of 100, 170, and 200 psf in 
R-134a heavy gas. Reynolds numbers based on 
model chord ranged from 2.4 to 6.5 million and 
these test, conditions corresponded to reduced fre- 
quencies (k) from 0.011 to 0.579 for the BSCW (1 
Hz to 30 Hz). Boundary-layer transition was fixed 
at 7.5 percent chord and the OTT pitch axis was 
located at x/c=0.3. The BSCW model has been 
previously tested at the TDT as part of the Bench- 
mark Models Program 4- ’ 43 during which a large 
database of unsteady pressures were obtained dur- 
ing motions on a flexible pitch and plunge mount. 
The purpose of this model was to use the data 
acquired using the OTT for the identification of 
experimental ROMs. 

Experimental Methodology 

Unsteady pressure measurements were made on 
the RSM and the BSCW while the models under- 
went pitch oscillations on the OTT at frequencies 
from 1 to 10 Hz (RSM) and from 1 to 30 Hz 
(BSCW). In addition, unsteady pressures were ac- 
quired during RSM/OTT and BSCW/OTT step 
inputs in order to provide data to compute aero- 
dynamic impulse responses. 

The identification of experimental unsteady 
aerodynamic (pressure) ROMs can be performed 
by using the same techniques used to identify the 
computational unsteady aerodynamic ROMs. The 
Volterra. theory of nonlinear systems is used a.s 
the basis for modeling the linear and nonlinear 
dynamic response of the unsteady aerodynamic 
system under investigation, as described in the ref- 
erences. 

For the present study, the identification of exper- 
imental unsteady aerodynamic impulse responses 
will be limited to the first-order, or linearized, 
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kernel. It is referred to as a linearized kernel 
since identification of the kernel (impulse response) 
may occur about a nonlinear steady-state condi- 
tion (such as St. transonic Mach number). Future 
research will focuiS: on the identification of the 
second-order kernel. The frequency-domain ver- 
sion of the second-order kernel is known as the 
bi-spectra, which has found important applications 
across a wide variety of disciplines for quantifying 
experimental nonlinear dynamics. 44 

The identification of the experimental unsteady 
aerodynamic impulse responses (first-order kernel) 
will consist of the deconvolution of a given in- 
put/output pair. Deconvolution involves the ex- 
traction of the impulse response of a system when 
the input and corresponding output are known. 
The input, in this case, is a sequence of positive 
and negative step inputs in pitch applied using the 
OTT and the output is any of several measured 
pressure responses from the wind-tunnel models. 
Deconvolution is then used to extract the impulse 
response for the given input/output pair. For the 
given OTT step input, an impulse response can be 
identified for each pressure measurement (sensor) 
on the wind-tunnel models. 

Once the impulse response has been generated, 
convolution is used to predict the pressure response 
due to sinusoidal inputs in pitch at various frequen- 
cies. 34 The measured results are compared to the 
predicted results (via convolution) in a subsequent 
section. 

Computational Results 

The AGARD 445.6 Aeroelast.ic Wing has been 
used extensively by several authors to validate 
computational methods. 13, 19,45 Although the 
aeroelast.ic behavior of this wing is fairly benign, 
the aeroelast.ic data, from the flutter test, of this 
wing has been used significantly in the past, for 
the validation of computational techniques. 45 The 
wing is a. 45-degree swept.-back wing with a. NACA 
65A004 airfoil section, panel aspect, ratio of 1.65, 
and a. taper ratio of 0.6576. The shapes of the 
first, four structural modes for this wing are pre- 
sented in Figure 6. The modes are first, bending, 
first, torsion, second bending and second torsion. 
The corresponding modal frequencies in vacuo are 
9.60, 38.2, 48.35, and 91.54 Hz. Additional details 
regarding this wing can be found in the references. 

AGARD 445.6 Wing 

Full CFD Flutter Solution 

This section presents results based on the tra- 
ditional full CFD flutter solution. The flutter 
solution is obtained by iterating between the non- 
linear aerodynamic system (flow solver) and the 
structural system at. a. given Mach number and dy- 
namic pressure entirely within the CFD code. The 
output, of this solution consists of a. time history 


Mode I 


Mode 3 Mode 4 

Fig. 6 Structural modes for the AGARD 445.6 
Wing. 

of the generalized coordinates of the aeroelast.ic 
system. Depending on the nature of this a.eroe- 
la.st.ic response (divergent, or convergent.), a. new 
dynamic pressure is selected and a. corresponding 
flutter solution is computed. This iterative process 
is used to define the flutter boundary at. several 
Mach numbers. The results presented in this paper 
are for the solution of the Euler equations within 
CFL3Dv6.0. 

Figure 7 presents the response of each of the 
four generalized coordinates at. a. Mach number of 
0.9, a. dynamic pressure of 89.3 psf, and a. struc- 
tural damping value (g) of zero. The divergent, 
nature of the first, and second modes indicates that, 
this condition is above the flutter boundary. By 
performing similar analyses at. different, dynamic 
pressures, a. dynamic pressure of 75 psf was deter- 
mined to be the flutter dynamic pressure (neutral 
stability point.) for this Mach number. The cor- 
responding flutter frequency was 14.8 Hz. The 
aeroelast.ic response at. a. dynamic pressure of 75 
psf is presented as Figure 8, indicating a. condition 
very close to the neutral stability of the a.eroe- 
la.st.ic system at. this condition. These solutions 
were computed using a. non-dimensional time step 
of 0.3 with 5 subit.era.t.ions per time step and us^ 
of mult.igrid capability for error reduction and con- 
vergence acceleration. Results from the full CFD 
flutter analysis 23 are consistent, with those from 
Lee-Ra.usch 13 and other Euler flutter results. 19 

The computational cost, for one flutter solution 
(at. a. given Mach number and dynamic pressure) 
is 71 CPU hours for the number of cycles shown 
in Figures 7 and 8. This is the total CPU cost, 
but., using 96 processors, the actual execution time 
is approximately 45 minutes on an Origin 2000 
cluster. The total time elapsed from the moment, 
the job is submitted for execution, however, can 
vary depending on the number of other jobs (from 
different, users) in the queue for the computer re- 
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Fig. 7 Aeroelastic transients in terms of gen- 
eralized coordinates at M=0.9 and q=89.3 psf 
for AGARD 445.6 Wing. 



Fig. 8 Aeroelastic transients in terms of gen- 
eralized coordinates at M=0.9 and Q=75.0 psf 
for AGARD 445.6 Wing. 

sources. The total elapsed time for a single flutter 
solution can therefore range from 45 minutes to 
several hours if the job queue is busy. In addition, 
although four cycles of the lowest frequency mode 
appear to be sufficient, for visually determining the 
stability of the system, accurate computation of 
the relevant aeroelastic frequency and damping re- 
quires additional cycles. If the number of cycles 
is doubled to eight cycles, the computational costs 
increase proportionately to 142 CPU hours and 90 
minutes of execution time. The total time elapsed 
can range from 90 minutes to several hours de- 
pending on the number of jobs in the queue. 

These costs are, of course, a function of the 
aeroelastic properties of the system under inves- 
tigation. The use of parallel processing clearly 
provides significant improvement in computational 
time. However, the computational costs (CPU) are 


still high because parallelization, obviously, does 
not reduce the amount of computation that needs 
to be done. Nonetheless, this is a significant im- 
provement over computations performed on a se- 
rial platform. 

Assuming four dynamic pressure solutions per 
Mach number, the cost of computing a flutter point 
(at one Mach number) is 568 CPU hours, requir- 
ing at least 360 minutes of execution time. The 
actual time invested, however, can be on the order 
of days since the value of dynamic pressure selected 
for the subsequent analysis depends on the results 
obtained from the previous analysis. If additional 
analyses involving parametric variations of struc- 
tural parameters (damping and frequencies) are 
needed, additional flutter solution^ would be re- 
quired, increasing computational costs (CPU and 
time). Finally, as can be seen, the output of tradi- 
tional CFD-based flutter analyses are aeroelastic 
transients which provide frequency and damping 
information at a given flight condition. These tran- 
sients can certainly be used to define the flutter 
boundary of the aeroelastic system under investi- 
gation but do not comprise a mathematical model 
of the system itself. In order to develop a math- 
ematical model of the system itself, a ROM is 
needed. 

ROM Flutter Solution 

Coupling the state-space model of the unsteady 
aerodynamic system with a state-space model of 
the structure within MATLAB/SIMULINK results 
in a state-space aeroelastic system. The aeroelastic 
response of the system is a function of the initial 
conditions of the structure and the dynamic pres- 
sure. 

In order to validate this state-space aeroelas- 
tic system, simulations were performed at vari- 
ous dynamic pressures. Figure ?? presents the 
generalized coordinate time histories and the cor- 
responding generalized coordinate FFTs at zero 
dynamic pressure (wind-off). The zero dynamic 
pressure eliminates all aerodynamic effects, leaving 
only structural effects. With zero structural damp- 
ing, the response consists of, in the time domain, 
the simple harmonic motion of the uncoupled vi- 
bration modes and, in the frequency domain, fre- 
quency peaks of the uncoupled vibration modes: 
9.60, 38.2, 48.35, and 91.54 Hz. 

At a dynamic pressure of 75 psf, Figure ??, flut- 
ter is evident. The corresponding flutter frequency 
of 14.5 Hz is clearly identified in Figure ??. A 
close-up of this aeroelastic transient is presented 
as Figure 11. This result shows excellent com- 
parison with the result of Figure 8, which was 
computed using CFL3Dv6.0 directly. In fact, the 
ROM results showed excellent comparison with 
results using CFL3Dv6.0 directly at all dynamic 
pressures investigated. 
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Fig. 9 Aeroelastic response of the reduced- 
order aeroelastic model at M=0.9 and q=0 psf 
for AGARD 445.6 Wing. 



Number of Time Steps 
a) Time domain 



b) Frequency domain 

Fig. 10 Aeroelastic transients in terms of 
generalized coordinates for the reduced-order 
model at M=0.9, q=75 psf, and g=0.0 for 
AGARD 445.6 Wing. 

These aeroelastic transients are computed in sec- 
onds within MATLAB/SIMULINlv, thus allowing 
a larger number of cycles to be computed for im- 
proved frequency resolution. In addition, if para- 
metric variations that involve the structure are 
desired (structural damping, updated frequencies, 
etc), the analyses can be performed since the un- 
steady aerodynamic system is unaffected by these 
variations. 

These results validate the ROM methodology 
presented and are examples of a new and pow- 
erful tool available to the aeroelasticia.n. Most 
importantly, the state-space models developed are 



Fig. 11 Close-up of the aeroelastic transients 
for the reduced-order model at M=0.9, q=75 
psf, and g=0.0 for AGARD 445.6 Wing. 

suitable for use within a mutidisciplinary design 
environment, including ASE analysis and design. 

Rigid Semispan Model 

Results comparing fully coupled aeroelastic so- 
lutions to ROM aeroelastic solutions using the 
softened version of the RSM' are presented Fig- 
ure 12, Figure 13, and Figure 14. Even though 
some differences in the responses is apparent, the 
ROM aeroelastic solutions are in excellent agree- 
ment with the fully coupled aeroelastic solutions. 
Both methods predict the neutral stability condi- 
tion (flutter) at approximately 123.8 psf. 

Each fully coupled GFL3Dv6.0 computation re- 
quired approximately 108 CPU hours on an Origin- 
3000 computer. Each impulse response (subse- 
quently used to develop aeroelastic ROM) required 
approximately 312 hours. The computational cost 
of the impulse responses could have been reduced 
by optimizing time step, subiterations, and other 
parameters. It is important to note, however, that 
each impulse response computation provides the 
wide range of reduced frequencies as well as the 
entire range of dynamic pressures. On the other 
hand, each fully coupled calculation only provides 
a single dynamic pressure. Research is currently 
focused on the simultaneous excitation of multiple 
modes in order to reduce the number of impulse 
responses required (currently one per mode) to 
generate the aeroelastic ROM. 

Experimental Results 

Rigid Semispan Model 

Results for the RSM/OTT are presented in this 
section. The RSM has a total of 131 in-situ pres- 
sure transducers. A step input in pitch using the 
OTT results in 131 unsteady pressure responses. 
Therefore, deconvolution can be applied to all of 
the unsteady pressure measurements to yield 131 
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c) Mode 3 

Fig. 12 Fully coupled versus ROM aeroelas- 
tic solutions for the softened RSM at M=0.7, 
q=118.08 psf. 


500 1000 1500 2000 2500 3000 3500 4000 
Time Steps 

c) Mode 3 

Fig. 13 Fully coupled versus ROM aeroelas- 
tic solution for the softened RSM at M=0.7, 
q=123.84 psf. 


unsteady aerodynamic impulse responses. For the 
sake of brevity and to demonstrate the feasibility 
of the method, results are presented for only one 
pressure measurement located on the upper sur- 
face of the RSM at the 60%) span location and the 
•30%) chord station. The data was acquired at a 
Mach number (M) of 0.8, a dynamic pressure (q) 
of 150 psf, and with the RSM at zero degrees angle 
of attack. 

Figure 15 presents the step pitch input com- 
manded to the OTT and the resultant pressure 
response at the pressure transducer location men- 
tioned above. Although a theoretical step input 
consists of an infinite slope where the step occurs, a 
physically realizable step input, such as that com- 
manded by the OTT, will be limited by the pitch 
inertia, stress, and load limitations of the model 
undergoing pitch. As can be seen, a step input that 
closely approaches a theoretical step input can, in 
fact, be applied by the OTT. The unsteady pres- 
sure measurement is also quite noisy, as can be 


seen. 

Using the sequence of step pitch motions of 
the OTT as the input and the unsteady pres- 
sure measurement as the output, deconvolution is 
applied to identify the unsteady aerodynamic im- 
pulse response. Figure 16 presents the time- and 
frequency-domain versions of the pressure impulse 
response identified via deconvolution. As can be 
seen in Figure 16(b), the identified impulse re- 
sponse exhibits significant frequency content, as is 
to be expected for an impulse response. An analy- 
sis of the unsteady aerodynamic impulse responses 
at all pressure transducer locations can provide a 
spatial mapping of the frequency characteristics 
of a given configuration at a given test condition. 
This type of spatial mapping may be useful for the 
design and optimal placement of various flow con- 
trol devices. 

Upon identification, the unsteady aerodynamic 
impulse response can then be used to predict the 
unsteady aerodynamic response due to any OTT 
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Fig. 14 Fully coupled versus ROM aeroelas- 
tic solution for the softened R.SM at M=0.7, 
q=129.60 psf. 
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a) Step Input in Pitch 



b) Pressure Response due to Step Input in Pitch 

Fig. 15 Commanded pitch motion and resul- 
tant pressure response on the upper surface of 
the RSM at 60% span and 30% chord at M=0.8, 
q=150 psf. 


input using convolution and the impulse response 
of Figure 16. In the following figures, comparisons 
are made between predicted unsteady aerodynamic 
responses arid the measured responses for several 
sinusoidal OTT motions. 

Figure 17 presents the comparison between the 
measured pressure response and the corresponding 
predicted pressure response for a commanded os- 
cillation of 1.2 Hz. The comparison is excellent and 
demonstrates the ability of the method to capture 
the dominant (driving) frequency while filtering 
out uncorrelated noise. The process of identifying 
the correlation for a given input/output pair, via 
deconvolution, has the added benefit that it filters 
out any information that is not correlated to the 
input /output correlated. Therefore, uncorrelated 
measurement noise, for example, is automatically 
removed as the impulse response is generated. This 
filtering capability is visible in Figure 17(b). 

Figure 18 presents the comparison between the 
measured pressure response and the correspond- 


ing predicted pressure response for a commanded 
oscillation of 10.1 Hz. For this case, without the 
predicted response, it would be very difficult, to dis- 
cern any periodicity in the measured response. The 
filtering capability of the deconvolution method 
proves to be essential at this frequency. 

At this condition, M=0.8 and q=150 psf, the 
linearity of the measured pressure responses (for 
this pressure transducer location) is defined by the 
excellent correlation between the experimental re- 
sults and the results computed using convolution. 
If predicted results do not compare well with mea- 
sured results, this could be an indication that some 
nonlinearity has influenced the measured response. 

In addition, because deconvolution involves in- 
put/output correlation, any uncorrelated white 
noise (measurement noise) is easily filtered out. 
Simple low-pass or high-pass filters would not be 
able to match this level of filtering capability and 
much more sophisticated band-pass filters would 
have to be introduced. However, even with band- 
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Fig. 16 Pressure impulse response obtained 
via deconvolution for the RSM; time domain 
and frequency domain (magnitude). 

pass filters, the question of which frequency range 
to filter would remain a serious question for the 
analyst,. Deconvolution automatically handles the 
filtering without a priori definition of a frequency 
range where filtering is desired. Analysis of 1 1 m m' 
results can subsequently be used to identify regions 
of linear and nonlinear behavior which will be help- 
ful in understanding dominant flow physics. 

Benchmark Supercritical Wing 

For the BSCW on the OTT, results are pre- 
sented for pressure measurements at the trailing 
edge, at a Mach number of 0.5, a dynamic pres- 
sure of 100 psf, and zero degrees angle of attack. 

Figure 19 presents the measured and predicted 
pressure responses due to a 2 Hz sinusoidal motion 
of the OTT in the time domain and frequency do- 
main (magnitude). The comparison is very good 
with excellent filtering of the uncorrelated noise. 

Finally, Figure 20 presents the measured and 
predicted pressure responses due to a 15 Hz si- 
nusoidal motion of the OTT in the time domain 



Time, secs 

a) Time domain 



10 1 10° io' io 2 


Frequency, Hz 

b) Frequency domain 

Fig. 17 Measured and predicted pressure re- 
sponses due to a 1.2 Hz sinusoidal motion of the 
OTT for the RSM; time domain and frequency 
domain (magnitude). 

and frequency domain (magnitude). For this case 
the comparison is good but Figure 20(a) indicates 
a difference in the response between the upper 
and lower portions of the measured response not 
captured by i lift, predicted response. Additional 
analyses are required to determine if this differ- 
ence is due to a possible nonlinear effect. Although 
at this subsonic Mach number nonlinear aerody- 
namic effects are not expected in general, since 
this pressure measurement is at the trailing edge, 
some local flow separation induced by the higher 
frequency may be occuring. 

Concluding Remarks 

The identification of computational and experi- 
mental reduced-order models (ROMs) for unsteady 
aerodynamic and aeroelast.ic analyses was pre- 
sented. Results included the implementation of 
the ROM technique into the CFL3Dv6.0 code, ere-* 
ation of a state-space model of the unsteady aero- 
dynamic system using CFD-ba.sed responses, and 
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Fig. 18 Measured and predicted pressure re- 
sponses due to a 10.1 Hz sinusoidal motion of 
the OTT; time domain and frequency domain 
(magnitude). 

integration of the state-space ROM into an aeroe- 
lastic simulation using the MATLAB/SIMULINK 
environment. The results presented for the iden- 
tification of an experimental ROM included re- 
sults for two wind-tunnel models: a Rigid Semis- 
pan Model (RSM) and a supercritical rectangular 
wing (BSCW). The experimental ROM identifica- 
tion method was described and used to predict the 
pressure responses due to various sinusoidal oscilla- 
tions of the Oscillating Turntable (OTT) at NASA 
Langley’s Transonic Dynamics Tunnel (TDT). The 
accuracy of the method was evaluated and the 
results presented can be used to gain insight re- 
garding the dominant unsteady flow physics. Ri 
addition, potential use of the method for data fil- 
tering applications was described. 
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